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Abstract. A systematic analysis of the A-dependence of 0-meson production in proton-nucleus collisions is pre- 
sented. We apply different formalisms for the evaluation of the 0-meson distortion in nuclei and discuss the theoretical 
uncertainties of the data analysis. The corresponding results are compared to theoretical predictions. We also discuss 
the interpretation of the extracted results with respect to different observables and provide relations between frequently 
used definitions. The perspectives of future experiments are evaluated and estimates based on our systematical study 



are given. 

PACS. 11.80.Fv Eikonal approximation - 
intermediate-energy reactions and scattering 



11.80.La Multiple scattering - 13.75.-n Hadron-induced low- and 
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1 Introduction 

The modification of hadron properties in a nuclear environ- 
ment remains one of the most mysterious problems in nuclear 
physics, see e.g. [1-7]. First experimental results on dilepton 
spectra from heavy ion collisions indicated a substantial mod- 
ification of the spectra in the vicinity of the p-meson mass [8- 
1 0] . Very recent results from PHENIX for high-energy Au+Au 
scattering [11] also indicated a significant enhancement of very 
low mass dielectrons within the energy range from 150 to 750 
MeV This modification of the high mass dilepton spectrum is a 
key question to understand the in-medium effects. For instance, 
the observation of a change of the 0-meson spectral function 
would unambiguously clarify the role of in-medium effects. 1 
Unfortunately, the statistical accuracy of the recent data [11] 
is still not good enough to draw definite conclusions about the 
modification of the 0-meson spectral function in heavy ion col- 
lisions. 

On the other hand there is a certain belief [12-14] that the 
</>-meson is almost not distorted in nuclear matter and thus can 
provide information about the early partonic stage of heavy ion 
collisions or the Quark Gluon Plasma (QGP). This means that 
0-meson production remains mostly unaffected by hadronic 
interactions and might serve as the perfect penetrating probe. 
Furthermore, it was proposed [12] that the formation of the 
QGP may be detected by enhanced 0-meson production result- 
ing from the absence of the Okubo-Zweing-Iizuka suppression 
that is quite substantial in elementary reactions [15-17]. In this 
case it is important that the distortion of the 0-meson is almost 



1 Note that the <ji>-meson is narrow and well isolated in mass from 
other mesons, which makes it a perfect probe of the in-medium dis- 
tortion. 



negligible. The most recent results obtained by the STAR Col- 
laboration [14] indicate some distortion of </>-meson production 
in Au+Au collisions, however, it is difficult to quantify the size 
of this effect. 

Some indications of ^-meson modification at normal nu- 
clear densities have been reported from measurements involv- 
ing photon [18] and proton beams [19-23]. The most remark- 
able observation in all these measurements is the anomalous 
A-dependence of </>-meson production. The modification of the 
spectral function implies the change of the pole position and 
the width. When the effective in-medium width of the 0-meson 
becomes larger, then its decay probability also increases, lead- 
ing to a stronger distortion in nuclear matter. Thus one might 
expect [24-29] that the A-dependence indicates a in-medium 
modification of the imaginary part of the </>-meson spectral func- 
tion. 

Our previous systematic analysis [29] of 0-meson photo- 
production from nuclear targets, however, showed that the ob- 
served A-dependence can be well understood by including the 
u)-</> mixing in nuclei. We have investigated single and coupled- 
channel phenomena and reproduced the recent SPring-8 data 
[18] as well as the older Cornell results at higher energies [30]. 

Here we analyze the available data on 0-meson production 
in proton-nucleus collisions. Most of the measurements were 
done at high energies. We provide the formalism for the data 
evaluation and give estimates for the theoretical uncertainties. 
Moreover, we review different formulas used in the analysis 
of the A-dependence at high energies and discuss their com- 
patibility. We also collect the different theoretical predictions 
relevant to the </>-meson distortion in nuclear matter. Finally, 
we hope that all this information allows for a systematic under- 
standing of the problem and helps in planning future strategies 
for experiment analysis. Therefore, we provide estimates for 



2 



A. Sibirtsev et al: ^4-dependence of 0-meson production mpA collisions. 



the data analysis currently carried out by the ANKE Collabo- 
ration at COSY [28]. 

2 Formalism 

Recently, the eikonal formalism has been considered to extract 
the in-medium properties of hadrons from the yl-dependence 
of their production on nuclei using photon and proton beams. 
The distortion of the produced hadrons at low nuclear matter 
densities can be expressed in terms of the so-called tp approx- 
imation. Here p is the nuclear density, while t stands for the 
forward scattering amplitude. This amplitude is frequently re- 
ferred to as /(0) and we follow this convention. This amplitude 
in nuclear matter is not necessary equal to that in free space be- 
cause it can be modified by the nuclear environment. Since this 
amplitude can be evaluated from the nuclear data, one might be 
able to extract the in-medium modification of the hadron inter- 
action. The tp approximation is well applicable at normal, but 
not high, nuclear densities and thus can provide a sensitive way 
for the nuclear data evaluation. See, e.g., Refs. [26,31-33] for 
evaluations using the tp approximation in the eikonal formal- 
ism. 

Originally this method was proposed to study the inter- 
actions of unstable hadrons with nucleons. Indeed, it is im- 
possible to produce a beam of unstable hadrons, like u, p, <fi, 
etc. mesons and measure their scattering on a hydrogen target. 
However, a nucleus can be used as a source of unstable hadrons 
as well as a target for their interaction and detection. In that 
sense, the eikonal formalism explores the nucleus as a labora- 
tory constructed at fm distances, i.e. some kind of "Fermilab". 

Furthermore, the eikonal formalism provides an unique op- 
portunity to evaluate the scattering length for the interaction of 
an unstable hadron with a nucleon. Indeed by measuring the 
A-dependence of the produced particle with different energies 
one can extract the scattering amplitude as a function of energy. 
Interpolating this amplitude to threshold, one obtains the scat- 
tering length. This scattering length could also be modified by 
the nuclear medium and its extrapolation to free space requires, 
unfortunately, model dependent assumptions. 

In many applications of the eikonal method the distortion 
of the produced hadron is factorized out of the full formalism 
in order to analyze the data. Furthermore, the interpretation of 
the evaluated results is given in terms of different variables. 
Thus it is difficult to compare different analyses and to link 
the simplified methods that were applied to the original one. 
For this reason, we briefly discuss the eikonal formalism for 
particle production processes and give the parameter set that 
can be extracted from the data analysis. 

2.1 Eikonal approximation 

The basic idea of the eikonal formalism is to express the in- 
teraction of a particle with a nucleus in terms of effective two- 
body interactions. To leading order, the total reaction ampli- 
tude can therefore be built up from a sum of amplitudes on a 
single nucleon. The corrections to the leading order come from 
multiple interactions. The eikonal approximation was proposed 
by Glauber [34,35] for coherent and incoherent scattering of 



hadrons from nuclei. The incident particle is assumed to inter- 
act independently with each target nucleon as it moves along a 
straight line trajectory through the nucleus. 

In the original formulation the incident and final particles 
are identical. Nevertheless, the formalism can be applied to 
production processes, i. e. when incident and final particles are 
different. The model was extended by Formanek and Trefil [36, 
37] for the case of resonance production in proton-nucleus col- 
lisions and generalized by Berman and Drell [38, 39], Ross and 
Stodolsky [40] and Drell and Trefil [41,42] for vector meson 
photoproduction on nuclei. The best known extension and ap- 
plication of multi-scattering theory to the coherent and inco- 
herent production of particles were given by Margolis [43,44] 
and Kolbig and Margolis [45]. 

In this section, we sketch the formalism of the eikonal ap- 
proximation necessary to derive the A-dependence of incoher- 
ent particle production off nuclei. Our aim is to illustrate how 
the elementary amplitude on a quasi-free nucleon is related to 
the total production amplitude on a nucleus. This will clarify 
to which extent one can extract the elementary amplitude from 
the experimental data. In the following, we will also show how 
one can interpret this amplitude with respect to different nu- 
clear phenomena. 

We consider a two-body interaction at impact vector b and 
describe the transition from the incident particle i to a final state 
particle / by the profile function f (b) defined as 

r</(b) = 2^k J A./V(q) e - 4qb , (1) 

where fa is the elementary i+N^f+X transition amplitude 
at momentum transfer q, while k is the projectile particle wave 
number. We also define profile functions Fa and Tjf similar to 
Eq. (1) with fa and /// denoting the amplitudes for i+N — ► 
i+X and f+N^f+X, respectively. 

Using Huygen's principle, the transition amplitude is given 
by the profile function as 

/i/(q) = £J d 2 br if (b) e^ h . (2) 

The distortion of particle i is calculated under the assumption 
that each target nucleon j fixed at transverse (sj) and longitudi- 
nal positions (zj) independently modifies the wave of particle i 
passing through the target nucleon by the factor 2 

i-r u (b- Sj ). (3) 

The distortion of the outgoing wave / is defined in a similar 
way. 

The nuclear profile function r^ib) describes the transition 
from an initial nuclear state \I) to the final nuclear state \F) as 

(F\r l A f (b)\i) = J2(F\ n [i-r l4 (b-s,)] 

I Zj<Zi 

m 

xr if (b- Sl ) [] [l-///(b-s m )] \I). (4) 

z m >zi 

2 This is an essential difference to cascade-like models that con- 
sider distortion on the basis of squared amplitudes. Our approach takes 
into account the interference between the distortion amplitudes and 
accounts for the quantum dynamics. 
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From this equation it is clear that the nuclear profile function 
accounts for the overall distortion of the incident particle i by 
nucleons in the target nucleus before the production of particle 
/ at the positions Zj<zi, for the transition i+N^f+X at zi, 
and for the distortion of the produced particle / at z m >zi. 

Eq. (4) indeed illustrates that the A-dependence of particle 
production is entirely given by the distortion of both the projec- 
tile and the produced particles. In principle, Eq. (4) summarizes 
the physical content of the eikonal approximation. 

The nuclear transition amplitude is now related to the 
nuclear profile function by Eq. (2). The cross section for pro- 
duction of the particle / in the i+A collision, summed over all 
nuclear final states \F), is finally given as 



d (jA _ I fA |2 
dt 2^\Jif\ ' 



(5) 



where t=—q 2 is the four- momentum transfer squared. 

The summation over the nuclear final states \F) can be 
evaluated using closure and approximating the many body tar- 
get wave function u j (r i , . . . , ) = 1 1) by the product of single- 
particle density functions as 



|u/(ri, ...,r A )\ 2 = Y[p A (ri). 

One can assume further that 

J d 2 s dz ra(b—s)p(s, z)~ J d 2 sru(b—s) 



(6) 



+ 00 

J dz p(b, z) 



-/«(o) 



2inT{h) 
kA ' 



(7) 



where fu(0) is forward i+N^i+N scattering amplitude and 
the optical thickness function T(b) is defined as 



+ OC 

T(h) = Aj dz PA {h,z), 



(8) 



with r 2 =b 2 +z 2 . Similar relations also hold for f+N — > f+N 
scattering. Furthermore, we use the relation 



l + /«(0) 



2i7rT(b) 
kA 



exp 



2iirf u (0)T(h) 
k 



(9) 



in order to evaluate Eq. (5). 

As a consequence, the yl-dependence of the / particle pro- 
duction in i+A collisions becomes a function of the single- 
particle density function p A and of the forward scattering am- 
plitudes for the i+N^i+N and f+N^f+N processes. The 
transition amplitude for i+N^f+N itself does not depend on 
A under the assumption that it does not depend on the Fermi 
motion, which in principle is different for different nuclei. This 
assumption is not satisfied for particle production off nuclei 
at energies below the reaction threshold in free space, where 
the Fermi motion is an essential part of the production process 
i+N^f+N. 



The imaginary part of the forward fa scattering amplitude 
is given by the optical theorem, which is a straightforward con- 
sequence of 5-matrix unitarity and 



k 

lm/ij(0) = — a,, 

47T 



(10) 



where oi is total cross section for the interaction of the particle 
i with a nucleon. Therefore, the forward scattering amplitude 
can be written as 



/ii(0) = — (i + ai)ai, 

47T 



(ID 



where ai=Re/(0)/Im/(0) stands for the ratio of the real to 
the imaginary part of the forward scattering amplitude. Here 
one should note that Oi and on in nuclear matter are not neces- 
sary the same as in a free space. Eq. (10) explicitly illustrates 
why the total Oi reaction cross section is used in the data evalu- 
ation. 3 A similar relation holds for the final particle scattering. 
Since the ratio a is a priori unknown, it is generally neglected 
and the differential cross section of Eq. (5) is finally given as 



da A 
dt 



da 



N 



dt 



-(N A + e), 



(12) 



where da N /dt is the elementary cross section for / particle 
production in a collision of the particle i with a nucleon, while 
N A is the effective number of target nucleons involved in the 
interaction, 



N A -- 



1 



d 2 b 



iT(b)_ 



-<r/T(b) 



(13) 



The first term of Eq. (12) is the leading order term, which 
describes the production preceded and followed by the distor- 
tion of initial and final particles in the nucleus. Since the pro- 
duction is considered on a single target nucleon i+N^f+X 
this process can be addressed as a direct production mecha- 
nism. If a i =af=cr then Eq. (13) can be written as 

N A = Jd 2 be- aT( ^T(h), (14) 

and for er=0 we have N A =A. So when initial and final parti- 
cles are not distorted by the nucleus, the A-dependence of the 
production process is a linear function of the atomic mass num- 
ber. In that case the nucleus is absolutely transparent, where the 
nuclear transparency is defined as 



T R = N A /A. 



(15) 



The second term denoted e in Eq. (12) contains the correc- 
tions due to the scattering of the i and / particles out of the 
forward direction and contributions from multiple interactions. 
The next-to-leading order terms need to be controlled in or- 
der to perform the theoretical evaluation of the data. For this 
reason, we provide here the results of numerical calculations 
which can not be found explicitly in previous studies. 



Note that Uj can be converted back to the imaginary part of the 
forward scattering amplitude or other variables related to the total re- 
action cross section, as will be shown below. 
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Table 1. The multiple scattering collision numbers N v calculated 
from Eq. (19) for a 12 C target as a function of v and total cross section 
a. The v=\ term is the effective number of target nucleons involved 
in the interaction at leading order, while the terms for v> \ come from 
the multiple scattering correction e of Eq. (12). 



a (mb) 


v=\ 


u=2 


u=3 


u=4 


v=5 


10 


8.12 


1.5 


.25 


.04 


.004 


20 


5.8 


1.8 


.57 


.16 


.038 


30 


4.3 


1.7 


.75 


.29 


.102 


40 


3.4 


1.5 


.79 


.39 


.176 


60 


2.3 


1.1 


.71 


.46 


.279 


80 


1.8 


0.8 


.57 


.41 


.301 


100 


1.4 


0.7 


.46 


.35 


.275 



The multiple scattering correction can easily be estimated 
under the assumption that the forward scattering amplitudes 
f(q) are the same for i+N^i+N and f+N^f+N elastic 
scattering. By parameterizing f(q) as 

f(q) = /(0)cxp(-B g 2 ), (16) 

with the slope parameter B, the correction e can be obtained as 



A 1 



1 



!/ = 2 



v H v exp[Bq*(l--)]N v , 



where 



1 + a 2 a 



(17) 



(18) 



16tt B 

The multiple scattering collision numbers N v are given by 

^ v = v\j d2b aVl e ~ ffT(b) [T(b)],y - (19) 

Furthermore, N\ equals the effective collision number Na 
given by Eq. (14) and represents the leading term of the mul- 
tiple collision series. Additional higher-order corrections given 
by e in Eq. (12) are determined by contributions N„ with v>\. 
Now /i can be estimated as the ratio of elastic to the total inter- 
action cross section and its value can not exceed one. 

Note that e reflects the systematic uncertainty in the appli- 
cation of the given formalism. To estimate e by Eq. (17) one 
needs to calculate N v . The N v terms calculated by Eq. (19) 
for carbon and lead targets are listed in Tables 1,2 for different 
cross sections a and v. The calculations were performed with 
a nuclear density function pa taken as a Wood-Saxon distribu- 
tion as 

PA(r) = l + exp[{r-R)/d\ > (20) 
using the density parameters [46] 

i?=1.28A 1/3 -0.76+0.8A- 1/3 fm, d=V3/7rfm. (21) 

The multiple scattering terms given in Tables 1,2 should be 
compared to the leading term N\=Na- It is clear that multiple 
scattering corrections to the A-dependence might be important 



when p is close to one. In principle, the ratio da A /da N could 
be different from the leading term N\=Na due to the e contri- 
bution. 

For instance, for pp scattering at high energies the elastic 
cross section does not exceed ~ 15% of the total cross section, 
while at COSY energies this number accounts for ~40%. This 
leads to large systematical uncertainties in the evaluation of the 
f+N^f+N forward scattering amplitude from the nuclear 
data collected at low energies. 

Furthermore, N v can not be considered as an effective num- 
ber of multi-nucleon clusters in nuclei and should not be con- 
sidered as an estimate for the production mechanisms involving 
the interaction of the particles with few nucleons. The multiple 
collision numbers N„ are attributed to elastic scattering of ini- 
tial (i) and final (/) particles before and after the production 
process i+N^f+N on a single nucleon. 

Finally, within an eikonal approximation, the A-dependence 
for direct particle production is given by total i+N and f+N 
cross sections and by the single density function qa- Since <7j 
and <jf both depend on the type of particle, i.e. photon, pion, 
nucleon, etc., as well as on its kinetic energy, the A-dependence 
also is a function of those degrees of freedom. The eikonal 
approximation does not include the A-dependence of particle 
production due to the ejectile emission angle. 4 

Now, the solid lines in Fig. 1 show the effective collision 
numbers Na calculated by Eq. (13) for carbon and lead nuclei 
as a function of <r,. The results are shown for c/=10, 40 and 
100 mb. The calculations indicate reasonable sensitivity to both 
(jj and Of. 

The arrows in Fig. 1 indicate the A 2 I 3 dependence that is 
frequently assumed for inelastic hadron-nucleus reactions. The 
A 2 ' 3 dependence is in general discussed in terms of the ab- 
sorptive interaction of the incident particle at the nuclear sur- 
face. However, this not an unique explanation as is illustrated 
by Fig. 1 . The A 2 ! 3 dependence might result from various com- 
binations of <7j and Of, which reflect quite different physics. 

The data analysis of the A-dependence of particle produc- 
tion from nuclei is frequently done in terms of an exponent a 



4 The main advantage the eikonal approximation offers is that the 
multidimensional equations reduce to a differential equation in a sin- 
gle variable. This reduction into a single variable is the result of the 
straight line approximation involved. 



Table 2. The multiple scattering collision numbers N v calculated us- 
ing Eq. (19) for a 207 Pb target as a function of v and total cross section 
a. The v=\ term is the first order effective collision number, while the 
v>\ terms are due to the multiple scattering correction e of Eq. (12). 



a (mb) 


v=\ 


v=2 


v=3 


u=A 


v=S 


v=6 


10 


59.9 


32.3 


15.0 


5.8 


1.9 


0.5 


20 


24.4 


17.8 


13.8 


9.6 


5.8 


3.1 


30 


13.5 


9.6 


8.5 


7.6 


6.2 


4.6 


40 


9.1 


5.9 


5.2 


4.9 


4.7 


4.3 


60 


5.6 


3.2 


2.5 


2.3 


2.3 


2.3 


80 


4.1 


2.2 


1.6 


1.4 


1.3 


1.3 


100 


3.2 


1.7 


1.2 


1.0 


0.9 


0.8 
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Fig. 1. The effective collision numbers iV,4 calculated by Eq. (13) for 
C and Pb nuclei and for different <r/ total cross sections as a function 
of Oi. The arrows indicate the A^-dependence. 



fitted to experimental results for a double differential cross sec 
tion as 

d 2 a A 



dndT 



cA c 



(22) 



where c is some constant. In general, the exponent a is eval- 
uated from data collected at different kinematical conditions, 
such as emission angles Q, kinetic energies T of produced par- 
ticles, incident beam energies and different kinds of projectile 
(i) and ejectile (/) particles. In the eikonal approximation, the 
variety of kinematical conditions can easily be classified by 
considering the A-dependence as a function of the total o; L and 
07 cross sections. However, one should not expect validity of 
the eikonal approximation at production angles away from the 
forward direction. 

Furthermore, the data analysis in terms of the A a function 
introduces additional systematical uncertainties, which can be 
well understood by inspecting Fig. 1. We consider the follow- 
ing example: For 07=10 mb we obtain a=2/3, as is indicated 
by the arrows in Fig. 1 after fitting both 12 C and 207 Pb data. 
This ^-dependence corresponds to o"i~39 mb for 12 C, as is 
shown by the arrow in the upper panel of Fig. 1 . At the same 
time the lower panel for a 207 Pb target shows that a=2/3 cor- 
responds to £Ti~22 mb. This large uncertainty in <ii is reflected 
in the standard deviation of the a slope and vice versa. To avoid 
this uncertainty in the theoretical analysis it is more useful to 
analyze the ratio of the production cross section measured with 
different targets. 

Nevertheless, for completeness we evaluate Eq. (13) as a 
function of both Oi and 07 for various nuclei and fit the calcu- 
lations by A a in order to determine the average slope a. This 
average slope a is shown in Table 3 as a function of ai and 
07. The calculations were done for C, Al, Cu, Ag, Au, and 
Pb nuclei. (We do not indicate the calculations for 0^=07=0 
that obviously result in a=\) The average slope parameters a 



Table 3. The average slope a of the A a -dependence fitted to the ef- 
fective collision numbers calculated from Eq. (13) for different cross 
sections <7i and 07. The calculations were done for C, Al, Cu, Ag, 
Au, and Pb nuclei. 





Oi (mb) 


(mb) 





10 


20 


30 


40 


50 


60 


70 


80 


10 


.85 


.66 


.56 


.51 


.48 


.46 


.45 


.44 


.43 


20 


.76 


.56 


.46 


.41 


.39 


.37 


.36 


.35 


.35 


30 


.71 


.51 


.41 


.37 


.35 


.33 


.33 


.32 


.32 


40 


.68 


.48 


.39 


.35 


.33 


.32 


.31 


.31 


.31 


50 


.66 


.46 


.37 


.34 


.32 


.31 


.31 


.30 


.30 


60 


.65 


.45 


.36 


.33 


.31 


.31 


.30 


.30 


.30 


70 


.63 


.44 


.36 


.32 


.31 


.30 


.30 


.30 


.30 


80 


.63 


.43 


.35 


.32 


.31 


.30 


.30 


.30 


.30 


90 


.62 


.42 


.35 


.32 


.31 


.30 


.30 


.30 


.30 



are afflicted with large uncertainties, e.g. a = 0.35 ± 0.18 for 
er,; = 40 mb and 07 = 30 mb. For this reason, we discourage 
fitting cross sections for a single nucleus and recommend to 
analyze cross section ratios instead. 

The eikonal approximation provides the ^"-dependence 
with l<a<0.3. As is shown in Table 3, the A 2 / 3 dependence 
can be observed under various conditions given by o^ and 07. 
For instance, an A 2 / 3 dependence is expected for <7t=0 and 
£>7=50 mb, which might correspond to the photoproduction of 
7r, p, lu and other mesons. Finally, the v4 2 / 3 -dependence can not 
be addressed as only due to the interaction of the incident par- 
ticle at the nuclear surface. It is also clear that an interpretation 
of the A 1 / 3 dependence could not be given in a unambiguous 
way. 



2.2 Angular dependence 

Since in the actual experiments the data are collected at some 
fixed angles or integrated over a certain angular interval it is 
important to estimate how much the A-dependence is effected 
by the ejectile angle. Such an estimate can be done using a 
quasi-classical approximation as is shown below. 

As a beam of particles i passes through the nucleus, its in- 
tensity is attenuated due to the scattering out of the beam direc- 
tion. Since particles can be removed from the beam because of 
both elastic and inelastic interactions with the target nucleons, 
the attenuation is determined by the distortion cross section Cj. 

The attenuation probability of an i particle passing through 
the nucleus at impact parameter b and longitudinal positions 
from —00 to z is then given by 



Si(h, z)= exp 



dz' pA(b, z') 



(23) 



and T z (b) can be considered as the linear nuclear density. 

Thus Eq. (23) is a semi-classical description of particle i 
attenuation in matter and one might argue that in this case 
tji should be taken as an inelastic or absorption cross section 
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rather than total reaction cross section. In that sense an emis- 
sion of a particle out of the beam trajectory can be also con- 
sidered a distortion. That is why er,; is discussed as a distortion 
cross section. 

The incident particle i interacts with target nucleon at a 
transverse b and longitudinal z coordinate and produces the 
final particle /. Let us to consider that / is moving along the 
line fixed at an azimuthal angle <j> and polar angle 9 with respect 
to the incident particle beam direction. The attenuation proba- 
bility of passing / through the nucleus in that case is given as 



s f (b,z,e, 



1 

Z7T 



-a A d£p(|r c |) 



(24) 



where the integration is performed along the path of the pro- 
duced particle / defined by 



Tt = (b + £ cos <j) sin 9) 2 + (£ sin c 



Of + (z + £cos60 2 . 

(25) 

Here we again assume that both initial and final particle move 
along the straight trajectories before and after the production 
vertex, but now the ejectile trajectory depends on / emission 
angles <f> and 9. 

Finally, the effective collision number can be evaluated by 
integration over the nuclear volume as 



N A - 



=jd 2 bdzp{b,z)Si(b,z)Sf{b,z,6,<f>), (26) 



and now depends on the production angles. It is easy to show 
that Eq. (26) reduces to the eikonal formalism given by Eq. (13) 
after an integration over the azimuthal angle cj) at polar angle 
9=0°. 




200 



200 
a, (mb) 

Fig. 2. a) The effective collision numbers calculated using Eq. (26) for 
C and Pb nuclei and <x, =40 mb as a function of 07. The solid lines 
show the results for emission angle 9=0°, while the dashed lines are 
the calculations for 0=150°. b) The ratio Npb/Nc as a function of 07 
calculated for angles 6=0°, 45° and 150°. 
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Fig. 3. The slope of the A a -dependence as a function of the effective 
cross section 07 shown for the different production angles: 9—0° (cir- 
cles), 9=45° (squares) and 0=180° (triangles). The results are given 
for the direct production mechanism. 



Now Fig. 2a) shows the effective collision numbers calcu- 
lated from Eq. (26) for C and Pb nuclei and (7^=40 mb as a 
function of 07. Here the solid lines indicate the results obtained 
for a final particle emission angle 9=0°, while the dashed lines 
show the calculations for 0=150°. 

The calculations indicate a quite strong angular dependence 
of Na- Fig. 2b) shows the ratio of the effective collision num- 
bers Npt,/Nc as a function of 07 for different angles: 9=0°, 
45° and 150°. Note that the ratios are almost the same for 
small emission angles O<0<45°. But the ratio becomes large 
at (9=150°. 

In addition Fig. 2b) illustrates the uncertainty in the eval- 
uation of <Tf from the ratios of data for / particle production 
from different nuclear targets. Namely, it is clear that for dif- 
ferent emission angles the Np^/Nc ratios roughly saturate at 
07 >40 mb. This actually means that model is insensitive to the 
value of 07 if it exceeds the limit of ~40 mb. 

For completeness, Fig. 3 illustrates the slope a of the A- 
dependence given by Eq. (22). The results are shown for the 
different production angles and as a function of 07. The cal- 
culations were done with <7i=40 mb. Note that a is saturated 
above cr/~40 mb, while the slope varies significantly at cr/< 
20 mb. The variation of a with the emission angle is almost 
negligible at #<45°. Because of the distortion of the incident 
particle the maximal value of a is below one. The minimal 
slope is close to ~0.3. 



3 A-dependence due to two-step 
production 

Fig. 3 shows that, neglecting the distortion of the final particle, 
i. e. for af =0, one might expect the maximal value for the slope 
of the A a -dependence around ~0.7, which is essentially driven 
by the distortion of the incident particle given by <Ji used in 
our calculations. Considering proton-nucleus interactions one 
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can use <7;=40 mb in the range of proton energies from ~3 to 
10 3 GeV. However, many experiments [19,47^-9,21,22] done 
with high energy proton beams indicate a slope a~l. This is 
the case for </>-meson production in pA collisions. 

This apparent discrepancy can be understood quantitatively 
in terms of multi-particle production at high energies. The pos- 
sible scenario is production of many pions that interact inside 
the nucleus and produce 0-mesons. Since the flux density of 
these pions or their multiplicity can be large and their ener- 
gies are above the TrN-^(f)N threshold, the probability of this 
process could be larger than the probability of the direct pro- 
duction considered previously. Indeed the pions are distributed 
through the whole nucleus and thus the ^-meson can be pro- 
duced over the full volume of the target. Therefore one can ex- 
pect that the A-dependence of the 0-meson production in such 
a case is proportional to A, while neglecting the final distortion. 

Quantitative estimates can be done for the two-step process. 
Assume that at some impact parameter b and longitudinal point 
z the incident energetic particle i produces some intermediate 
state j. Due to the Lorentz boost, the j particle is moving along 
the beam direction, i. e. at the same impact parameter b and 
at some point z produces the final particle /. The final particle 
is now moving along the path given by the polar angle 9 and 
azimuthal angle <f> of the emission. 

Thus we have two sub-processes and this is called a two- 
step mechanism. The first one is iN^jN and the second one 
is jN^fN where the distortion of the i, j and / particles is 
taken into account. The probability of the first process can be 
derived in analogy to Eq. (23) and is given as 



z 

Sij(b,z) = \f l] \ 2 J d~z p A (b,z) 

— OG 

Z Z 

-Gi J dz f pA(h, z') - GjJ dz f paO>, z' 



exp 



, (27) 



where is the amplitude of the iN^jN transition, while 
(Ti and a-j account for the distortion of the i and j particle. 
The attenuation probability of passing / through the nucleus is 
similar to Eq. (24). 

The effective collision number for the two-step process is 



N A 



= J d 2 bdzp A (b,z)S ij (b,z)Sf(b,z,6,<f>), (28) 



where Sf is given by Eq. (24). Note that for Sij=Sf=l the 
A-dependence is proportional to A. 

Now Fig. 4 shows the slope a of the A a -dependence as 
a function of the effective cross section a/. The calculations 
were done using Eq. (28) and the results are shown for the 
different production angles. Note that the slope for the two- 
step production mechanism can be substantially larger than for 
the direct one. For the two-step process a can be close to one. 
However one can not estimate the /4-dependence due to the 
multi-step production and thus the extraction of <jf can not be 
done unambiguously. 

A quite different estimate of the ^-dependence was pro- 
posed in Ref. [50] for the evaluation of charmonium absorp- 
tion from high energy J/if'-production in p+A collisions. In 
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Fig. 4. The slope of the A a -dependence as a function of the effective 
cross section 07 shown for the different production angles: 9—0° (cir- 
cles), #=45° (squares) and #=180° (triangles). The results are shown 
for the two step production mechanism and the calculations were done 
using Eq. (28). The solid line indicate the results obtained by Eq. (29). 
The dashed line is the (pL) approximation calculated for Ag nucleus. 
The shaded box indicates the experimental result from the HERA-B 
Collaboration. 



that case the ^-dependence is given by 



00 

N A = Jd 2 bdz pA{b, z) exp ajj dz' PaO 3 , z 1 ) 



•(29) 



Following this formalism it is clear that if there is no distortion 
of the final particle the A-dependence is proportional to A and 
this is the limit, which might be expected from the nuclear data 
on incoherent particle production. The solid line in Fig. 4 indi- 
cates the slope a of the A a dependence calculated by Eq. (29) 
for different <r/. Note that in the derivation of Eq. (29) it was 
assumed [50] that (if is small, so one can not seriously discuss 
the difference between the two-step formalism and the above 
estimates at large distortion cross sections. 

Another simple evaluation of the A-dependence can be done 
using so-called pL parameterization [51] 



N A = A exp [~a f (pL)} 



(30) 



where (pL) is the average amount of matter crossed by the final 
particle and 5 



(pL) 



A-l 
2A 2 



d 2 b[T{h)} 2 



(31) 



Now if a/ is small one can use the following expression 

(PL) 



N A = A a , a = 1-07 



In A 



(32) 



Note that our normalization of the nuclear density is A, while in 
Refs. [50,51] it is unity. 
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which was extensively applied in the evaluation of the distor- 
tion of the charmonium cross section in nuclear matter. 

Apparently the slope a can be calculated using only one tar- 
get and as we found only slightly depends on A unless one uses 
light targets. The dashed line in Fig. 4 shows the a obtained by 
this (pL) approximation using Eq. (32) and an Ag target. As 
was mentioned before, the approximation is valid for small a f 
and indeed is in rough agreement with Eq. (29) for <7f<20 mb. 

To obtain experimental values for the slope a one needs the 
production cross sections measured for different nuclear targets 
A. It was found [51] in the evaluation of charmonium absorp- 
tion, that the value of a extracted from a fit to a given data set 
depends on the nucleus used as the lightest target. Indeed the 
experiments that use heavy targets with hydrogen or deuterium 
systematically obtain large values of the slope a. 

4 0-meson production in p+A collisions at 
high energies 



Table 4. The slope a of the ^"-dependence of </>-meson production 
obtained in the different experiments and the distortion cross section 
(70 evaluated by Eq. (29). For the KEK-PS E325 results the first error 
in a is statistical and the second error is systematic. The results are 
shown for the different ranges of rapidity y and transverse momentum 
Pt given in GeV/c. 



Experiment 


Ref. 


a 


o> (mb) 


HERA-B 


[19] 


0.96±0.02 


2.1±1.2 




BIS-2 


[20] 


0.81 ±0.06 


12±4 


ACCMOR 


[21] 


0.96±0.04 


2.1 ±2 




NA 11 


[22] 


0.86±0.02 


9±2 


KEK-PS E325 


[23] 






y 


0.9-1.1 




0.916±0.101±0.022 


4.9±4 


y 


1.1-1.3 




1.050±0.101±0.02 


0±2.8 


y 


1.3-1.5 




0.881±0.084±0.02 


7.2±5.8 


y 


1.5-1.7 




0.780±0.119±0.019 


14±8.3 


pt 


0-0.25 




0.971±0.101±0.019 


1.7±7 


pt 


0.25-0.50 




0.890±0.066±0.019 


6.7±4.9 


pt 


0.50-0.75 




0.924±0.111±0.021 


4.4±4 



As was mentioned previously the results on A-dependence of 
inclusive (/.-meson production from p+A collisions at high beam 
energies indicate a large slope a~ 1 . Here we shortly review the 

current status. Moreover, we evaluate the (/.-meson distortion the dlfferent laboratory momenta of the produced (/.-mesons, 

cross section and collect the results in Table 4. Furthermore the That would 8 lve information about the momentum dependence 

interpretation of a, is given in the next Section. of the distortion. In spite of that uncertainty in the analysis of 

Most recently the ^-dependence of inclusive ./.-meson pro- the KEK " PS E325 data our ™ ult * for the </>-meson distortion 

duction off nuclei using a 920 Ge V proton beam was measured cross sectlon are summarized in the Table 4. 
with HERA-B detector at HERA storage ring [19]. This exper- The dependence of the inclusive 0-meson production in 

iment was done with C, Ti and W targets and the cf)—>K + K~ neutron-nucleus interactions at 30-70 GeV was studied by the 

decay mode was used for the (/.-meson reconstruction. The data BIS " 2 Collaboration at the Serpukhov accelerator [20] . Here 

analysis shows the slope a=0.96±0.02. As discussed before the C ' Al and Cu tar 8 ets were used and 11 was found that the 

only the multi-step mechanism can be an explanation of the slo P e a = 0.81±0.06. That corresponds to a distortion cross 

HERA-B observation. Indeed the shaded box in Fig.4 indicates sectlon of 1 2±4 mb. 

the result from the HERA-B Collaboration, which in princi- In Ref - i 21 l the ^-dependence of </> meson production by 

pie can be explained assuming cr =2. 1±1 .2 mb. Here we use a 100 GeV/c P roton beam was determined through the analysis 

Eq. (29) as was done in the analysis of J/& distortion. of the data collected with H 2 and Be targets. The measure- 

A systematic study of (/.-meson production in p+A colli- ments was done by the ACCMOR Collaboration at SPS. The 

sions at a beam energy of 12 GeV was carried ot by the KEK- <t>^K+K decay mode was used for the reconstruction. It was 

PS E325 Collaboration [52,23,53,54]. The slope a of the A a - found that slo P e a=0.96±0.04. Moreover, it was argued [21] 

dependence was evaluated using C and Cu targets. The most that the use of the #2 target in general introduces additional 

recent results [54] allows to investigate how a depends on the systematic uncertainties, which are difficult to estimate. This 

(/.-meson momentum as well as to obtain a for the (j)—>K + K~ result is close t0 me HERA-B observation, 
and 0->e+e _ decay mode. It was found that a is statistically The ^-dependence of the inclusive (/.-meson production from 

the same for these two different decay modes in the same kine- beryllium and tantalum targets using a 120 GeV proton beam 

matical region. Furthermore it turns out that a depends on reac- was studied with the NA 1 1 spectrometer at CERN SPS [22] . 

tion kinematics. Here we would like to make some comments. The data analysis indicates that a=0.86±0.02. Applying Eq. (29) 

The relevant kinematics for the evaluation of A-dependence one can estimate cr^~9±2 mb. 
is given by the final particle production angle. It is clear that Below we also list experiments that did not measure the 

this angle defines the path of the particle and therefore the A-dependence, but assume some values of a under certain as- 

amount of matter involved in the distortion. Eq. (29) is applica- sumptions in order to analyze the data. 

ble at forward angles, while Eq. (28) can be used for large an- Inclusive (/.-meson production off beryllium nuclei by 70 

gles but accounts only for the two-step production mechanism. GeV/c protons was studied with the Sigma spectrometer [47] 

Nevertheless within such limitations one can realize from Fig. 4 at the Serpukhov accelerator. For the evaluation of elementary 

that the angular dependence is essential for data evaluation. pN^(f>X cross section it was assumed that the A-dependence 

Unfortunately, the KEK-PS E325 data are given either as of the nuclear cross section is proportional to A - 7 , as was mea- 

a function of 0-meson momentum alone [54] or as a function sured for K ~ A interactions [48]. 

of rapidity and transverse momentum [23]. In our opinion, the High statistics (/.-meson production from p+Be collisions 

ideal case is to fix forward angles and to extract the slope a for at beam momenta of 120 and 200 GeV/c was studied by the 
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ACCMOR Collaboration at SPS [49]. The data evaluation was 
done under the assumption that the A-dependence is propor- 
tional to A, which was motivated by the experimental results 
published in Refs. [21,22]. 

Finally the evaluated distortion cross sections are collected 
in Table 4. Unfortunately large uncertainties in the experimen- 
tal results for a produce large uncertainties in a^. In our opin- 
ion, the analysis of the ratios of production cross sections 
from different nuclei with respect to the C-target results might 
be less uncertain. In that case the systematical errors might can- 
cel up to large extent. However, such an analysis requires mea- 
surements with many different nuclear targets, which is not the 
case for some experiments available now. 



5 Interpretation of <r f 
5.1 Definitions 

The interpretation of oi is a general problem. Following our 
derivation given in Section 2, the eikonal formalism operates 
with the forward scattering amplitude /(0) and 07 appears 
through the optical theorem. In that sense aj is the total cross 
section for the interaction of a particle / with a nucleon em- 
bedded in nuclear matter. In the classical derivation 0$ is con- 
sidered a distortion cross section. There is no conflict between 
these two definitions if we consider attenuation of the flux of 
the final particle due to all possible processes available in the 
nucleus. That might be general absorption, scattering out of the 
initial trajectory, decay of an unstable particle followed by the 
distortion of the decay products, interaction with few nucleon 
configurations and whatever one can assume. The total sum 
over all these possible processes is an effective total or distor- 
tion cross section. 

As we emphasized previously o$ is not the free vacuum 
f+N total cross section since it has to be extracted from the 
nuclear data and might be modified by in-medium effects. How- 
ever it is always worthwhile to compare 07 with the free cross 
section if that is available. That comparison would show whether 
additional reaction channels were open in nuclear matter or 
whether reaction channels available in free space are blocked 
in the nucleus. For instance some transitions might be blocked 
due to the Pauli principle. In the analysis of charmonium prop- 
erties in nuclear matter, the distortion cross section <jf is a stan- 
dard variable generally used everywhere throughout the rele- 
vant discussions. 

Since in some calculations not the distortion cross section 
but other variables are used we provide here some useful rela- 
tions for the conversion. Let us first remind the reader that the 
complex forward scattering amplitude is related to the cross 
section by Eq. (1 1). At the same time the complex local poten- 
tial is given in terms of the complex forward scattering ampli- 
tude /(0) as [55-58] 



called in-medium collisional width Ar and mass shift Am of 
the / particle, which are [59-62] 



V 



mjvrn/ 



(33) 



where p is local nuclear density. The potential depends on / 
due to the energy dependence of /(0). It is possible to use a so 



(34) 



Am = -2n mN + mf pRcf(0), (35) 

where to a? is the nucleon mass and kf is the momentum of 
the final particle. 6 In principle one can replace masses by total 
energies and discuss the Ar and Am at high energy of the 
final particle. Note that the in-medium collisional width and 
mass shift are not invariants and can be changed by a Lorentz 
boost, so one should use these variables in the rest frame of the 
/ particle. 



5.2 Estimates for er^ in vacuum 

It is useful to compare distortion cross sections extracted from 
the nuclear data with its values in free space. That allows to 
inspect directly the possible in-medium modification of the a^. 
There are various well known methods to estimate the <p+N 
interaction cross section. 

The <fi+N cross section can be evaluated in the Vector Dom- 
inance Model from the jN^(f>N reaction. Within VDM the 
hadron-like photon [63] is a superposition of all possible vector 
meson states. Therefore the jN^<j)N reaction can be decom- 
posed into the transition of the photon to a virtual vector meson 
V followed by the elastic or inelastic vector meson scattering 
on the target nucleon and production of the final 0-meson. The 
reaction amplitude is then written as [64,65] 



JjN^cfiN — JVN^<pN 



(36) 



where the summation is performed over vector meson states. 
Moreover, a is the fine structure constant, jy is the photon 
coupling to the vector meson V and fvN^<pN is the amplitude 
for the VN^xfiN transition. 

The coupling -fy is given by vector meson decay into a 
lepton pair [66] 



Amf 



1 



2mf 
ml 



(37) 



where my and mi are the masses of vector meson and lepton, 
respectively. Taking the di-electron decay widths [67], the pho- 
ton couplings to the lightest vector mesons are 



7p =2.51, 7^=8.47, 7^=6.69. 



(38) 



Note that non-diagonal, i.e. pN^<f>N and ujN^4>N as well 
as diagonal cpN-^cpN transitions contribute to the reaction am- 
plitude of Eq. (36). VDM suggests that the virtual vector me- 
son stemming from the photon becomes real through the four- 



6 It is clear that the introduction of Ar in the calculation described 
above requires an accurate definition of the nuclear density depen- 
dence rather then average estimate of Eq. (34). 
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Fig. 5. The <f)N cross section as a function of </>-meson total energy. 
The open circles show the results obtained by Eqs. (39, 40) from the 
data on forward 0-meson photoproduction cross section assuming s- 
wave dominance in the <f>N— ><j)N scattering. The closed circles are 
the estimates for cr^jv x yfl + a| given by Eq. (42). The squares are 
the results extracted [29] from the data on ^-photoproduction from nu- 
clei [30, 18]. The solid line is constrained by the additive quark model 
of Eq. (44). The dashed line shows the estimate using Eq. (43). 



momentum t transferred to the nucleon, which in general re- 
quires the introduction of a form-factor at the interaction ver- 
tices [68-70]. In many analyses [64], this form factor is ne- 
glected. Thus the VDM analysis of photoproduction data re- 
quires additional assumptions. 

There are many precise data on </>-meson photoproduction 
differential cross sections at energies close to the reaction thresh- 
old. These data can be used for evaluation of the (f>N^cj>N 
scattering amplitude squared applying the VDM as 



da 7 N^(pN 7T a 2 

2 2 \}<t>N^<t>N\ 



dt 



(39) 



Now if the scattering is dominated by s-waves one could ex- 
tract the <f)N scattering length at threshold. It is related to the 



cross section as 



&4>N — 47T \f^N^4>NY 



(40) 



and is shown in Fig. 5 by open circles. Here we use the for- 
ward 0-meson photoproduction cross section available near the 
reaction threshold [71,72]. This scattering length can be com- 
pared with other theoretical predictions. For instance the esti- 
mate based on QCD sum rules [73] provides a real tpN scatter- 
ing length of a0N~-O.15 fm that corresponds to a cross section 
of <70~2.8 mb and seems to be in good agreement with the data 
evaluated by Eqs. (39, 40) under assumption of s-wave <j)N 
scattering. However, note that the 0-meson photoproduction 
differential cross sections are essentially anisotropic already at 
energies close to the threshold and can be well parametrized 



7 In our normalization the scattering amplitude / equals to the scat- 
tering length a^jv at q<j,— »0. 



as da I dtoc exp (bt) with a slope b~3 GeV 2 and t being the 
four momentum transfer squared [71,72]. In that sense, the es- 
timates shown by an open circles in Fig. 5 might not be correct 
and should be taken with a grain of salt in the evaluation of <fiN 
scattering length. 

The estimate [74] based on an effective Lagrangian ap- 
proach predicts the scattering length a0Ar = (— 0.01+i0.08) fm 
that corresponds to a^O.8 mb. 

Another estimate can be obtained from a QCD van der 
Waals potential calculation [75], which predicts a </>-nucleon 
bound state. In that case the Born approximate scattering length 
is given by the potential as [76,77] 



a<f>N 



m N + ?Ti0 



dr r 2 V(r), 



(41) 



where the potential was taken in the Yukawa form V(r) = 
— aexp[— r/j]/r with strength a=1.25 and range /i=0.6 GeV. 
These parameters were obtained for an attractive potential and 
result in a0jv~O.67 fm. They correspond to 0^=56 mb at q^O. 

Furthermore, applying Eq. (11) the -fN^4>N differential 
cross section of Eq .(39) at t—0 can be written as 



da 



'fN-Kj>N 



dt 



4 



<pJ u (pN- 



(42) 



Since the ratio of the a$ and <f>N total cross sections are 
unknown, one can extract from the photoproduction data only 

\, which is shown by 



their combination, i. e. a^N x y 1 + 

closed circles in the Fig. 5. Here we use the data collected in 
Ref. [29]. If one assumes that a$=0 these results may be con- 
sidered as the energy dependence of the (f> N cross section. While 
for many processes the real part of the scattering amplitude 
vanishes at high energies, the ratio of real to imaginary part of 
the amplitude a is large at low energies and moreover substan- 
tially depends on the momentum of the scattered particle [67]. 
Just to illustrate such a possibility the dashed line in Fig. 5 
shows the dependence 



a (f)N x ^Jl + a 2 = 10 (mb) x 



0.6 (GeV/c) 



(43) 



Again this might ensure that at high energies the (j)N cross sec- 
tion approaches some value around 10 mb, but still does not 
provide a trustworthy estimate of a^N close to threshold. At 
least it is not appropriate to estimate the real part of the for- 
ward scattering amplitude and to evaluate the in-medium mass 
shift using Eq. (35). 

Within an additive quark model the (j)N cross section is 
given as [78] 



0>AT — &K-N + &K+N — OV-JV) 



(44) 



where the elementary cross sections are taken at the same in- 
variant collision energies. Since the 4>N reaction threshold is 
mA+mjv— E96 GeV one can safely use Regge parametriza- 
tion for the continuum or non-resonant meson-nucleon scatter- 
ing amplitudes [79]. Now Eq. (44) is shown by the solid line in 
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Fig. 5 and is in reasonable agreement with VDM results at high 
energies. 

It is also worthwhile to show the estimate based on the dy- 
namical study of the <j>N bound state within the chiral SU(3) 
quark model. By solving a resonating group method based equa- 
tion [80] it was found that the binding energy of the state might 
range from 1 to 9 MeV. With respect to the s-wave scattering 
length a^N, the relation between the pole of the S'-matrix and 
binding energy e is given as 



a>4>N 



2m,N rn<i> 

TCI TV + 7710 



-1/2 



(45) 



so that the real part of the scattering length ranges from 2.1 to 
6.3 fm. This scattering length is large compared to the other 
results. Note that this result is used in three-body calculations 
of the <f)NN nuclear cluster binding energy [81]. 



5.3 Estimates for <t in matter 

Only some results evaluated from high energy proton-nucleus 
collisions (summarized in Table 4) are in agreement with the 
data shown in Fig. 5. The uncertainties of the KEK-PS data [23] 
are still too large to draw a definite conclusion right now. The 
results from BIS-2 [20] and NA-11 [22] are consistent with 
vacuum estimates at high energies. 

However the results fromHERA-B [19] and ACCMOR [21] 
indicate a substantially smaller 0-meson distortion cross sec- 
tion. This observation is difficult to interpret since in high en- 
ergy experiments, the ^-mesons are produced with high mo- 
menta and should be almost blind to any in-medium modifica- 
tion. 

The squares in Fig. 5 show the 0-meson distortion cross 
section extracted in [29] from the data on </>-photoproduction 
from nuclei [30,18]. As we already discussed, it is not nec- 
essary that these in-medium results are the same as ct^n in 
vacuum. However we observe reasonable agreement between 
nuclear results and those evaluated by VDM at high energies. 

Substantial modification of slow 0-mesons in nuclear mat- 
ter was proposed in Ref.[74]. It was found that the mass of the 
</>-meson almost does not change in matter, while the change 
of the width accounts for Z\_T~45 MeV at normal nuclear den- 
sity. Following Eq. (34) one can estimate the distortion cross 
section as O0~7O mb for fc0=lOO MeV/c. 

The energy dependence of the in-medium 0-meson width 
was studied in Ref. [24]. While for E^—m^ the width is about 
20 MeV at normal nuclear density, it increases up to 40 MeV 
at a 0-meson energy of 1 . 1 GeV So it is really changed by a 
factor of two over 80 MeV in energy. 8 This corresponds to a 
variation of the distortion cross section from ~27 to 15 mb and 
seems to be in agreement with the dashed line shown in the 
Fig. 5. 

Furthermore, the energy dependence of the </)-meson width 
at normal nuclear density was investigated in Ref. [25]. It was 
found that the in-medium width Z\_T~22- 1 7 MeV slightly varies 
with energy within the range E^—ra^ to 1.2 GeV. 




Here we refer to the energy dependence of the <^-meson width in 
nuclear matter at different densities, which is shown in the Fig.4 of 
Ref. [24]. 



200 
(grad) 

Fig. 6. The ratio of the effective collision numbers calculated for 64 Cu 
and 12 C nuclear targets as a function of ^-meson production angle 
shown for the different within the range from to 40 mb with a 
step size of 5 mb. 



6 Predictions for </>-meson production in 
P+A collisions at COSY energies 



The COoler SYnchrotron (COSY) at Jiilich provides an unique 
opportunity to study the 0-meson distortion in nuclear matter 
at low energies. Our analysis indicates that even the vacuum 
(f>N interaction is not well understood and different estimates 
illustrated in Fig. 5 are in substantial disagreement at 0-meson 
energies below 2 GeV. Moreover, the available predictions [74, 
24,25] state that the in-medium modification of the 0-meson 
width is substantial at low energies, although the real size of 
that change is not well established. 

Such a situation requires precise measurements of 0-meson 
production from p+A collisions at low energies, as was pro- 
posed in Refs. [26,27]. A dedicated experiment on 0-meson 
production from the proton interaction with 12 C, 108 Ag and 
197 Au targets at maximum COSY energies of 2.83 GeV was 
proposed by ANKE Collaboration [28]. Here we show the re- 
sults for the A-dependence of </>-meson production in pA col- 
lisions at few GeV energies. Note that at low energies the <fi- 
meson production due to multiple processes is suppressed due 
to the final particle multiplicities and large 0-meson production 
threshold. 

For the further calculations we fix (7^=40 mb, which stands 
for the average cross section for the interaction of the beam 
protons with target proton and neutron. Although the pN in- 
teraction can be modified in nuclear matter one would not ex- 
pect that this effect is significant for the protons with momenta 
above ~1 GeV/c. 

We believe that the analysis of the ratios R of the produced 
0-meson contains less theoretical uncertainties than the analy- 
sis of the differential cross section d 2 a/dT /dfl itself for each 
nuclear target A or the slope a, as was discussed previously. 
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Thus in the following we show our predictions for the ratio 



R = 



d 2 a A 
dTd cos 



d 2 a c 



dTd cos 9 



Na 
N c 



(46) 



taken with respect to the carbon target. Here T is the kinetic 
energy and 8 is the emission angle of the produced 0-meson, 
while Na is the effective collision number that was calculated 
for the different and 6. Moreover, the analysis of the ratios 
has additional advantages since systematic experimental uncer- 
tainties can be substantially reduced. 

Once more we would like to emphasize that in the evalua- 
tion of an effective collision number Na we use an effective in- 
medium cross section cr/=cr^. This is not the cross section for 
the 0-meson interaction with free nucleon. Moreover, as was 
discussed previously the multiple scattering series corrections 
e given by Eq. ( 1 2) can not be isolated and thus the extracted o§ 
cross section contains such a multiple scattering contribution. 
Nevertheless it is of great importance to compare evaluated 
from the nuclear data with the vacuum 4>N cross section. 

The calculations were done for the 12 C, Cu, Ag and 
196 Au nuclear targets, which is in line with the targets pro- 
posed [28] for the measurements at COSY. Figs. 6-8 shows the 
calculated ratios as a function of the 0-meson production an- 
gle. The lines indicate the results for different <j<p given within 
the range from to 40 mb with a step size of 5 mb. 

The calculations indicate substantial angular dependence 
of the ratio. As we showed previously, the analysis of the A- 
dependence for the forward particle production is most prefer- 
able for several reasons. At forward angles, the reaction mech- 
anism can be formulated within an eikonal basis and contains 
less theoretical uncertainties than the quasi-classical approxi- 
mation. Furthermore, due to the scattering dynamics any pos- 
sible multiple processes contribute less to the forward particle 
production. 



o 




Fig. 7. The ratio of the effective collision numbers calculated for 
10S Ag and 12 C targets as a function of the 0-meson production an- 
gle shown for different a$ within the range from to 40 mb with a 
step size of 5 mb. 
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Fig. 8. The ratio of the effective collision numbers calculated for 
196 Au and 12 C targets as a function of </>-meson production angle 
shown for the different within the range from to 40 mb with a 
step size of 5 mb. 



Moreover, at forward angles the ratios indicate reasonable 
sensitivity to the distortion cross section when er^ stands below 



~20 MeV. At large angles and for large < 
very high precision data. 



, the analysis requires 



For completeness let us to illustrate how to use the figures 
with the calculated ratios. We take as an example the ratio of the 
effective collision numbers from Ag and 12 C nuclei shown 
in Fig. 7. Let us consider the production at forward angles, i. 
e. 6<10°. If there is no distortion of the incident proton and 
final 0-meson the ratio equals that given by the target mass 
numbers leading to R=9. Due to the distortion of the incident 
proton and neglecting the distortion of the produced 0-meson 
one finds that i?=4.6. That is the maximum value given by the 
direct 0-meson production mechanism. The A 2 / 3 dependence 
corresponds to R=43 and leads to a 0-meson distortion cross 
section of less than 5 mb. The A 1 / 3 dependence results in a 
ratio of ~2.1 and corresponds to >40 mb. 

Now the question arises if the measured ratio is larger than 
R=4.6. This explicitly indicates the contribution from multiple 
or two-step processes, as is illustrated by Fig. 4. In that case the 
extraction of the 0-meson distortion in nuclear matter is much 
more model-dependent. Then one could use additional kine- 
matical constraints in order to isolate direct production mecha- 
nism, as was done for instance in Ref. [82]. 

Finally, we can estimate the A-dependence in case the dis- 
tortion cross section is (7^=10 mb. Then one might expect at 
forward angles the ratio for ^-meson production from 108 Ag 
and 12 C targets to be R=3.2 as is illustrated by Fig. 7. This 
corresponds to a mass dependence of ocA - 52 . This result is 
compatible with in-medium width of ~30 MeV for an average 
0-meson momentum of 500 MeV/c. 



A. Sibirtsev et al: ^4-dependence of 0-meson production mpA collisions. 
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7 Conclusions 

A systematic analysis of the yl-dependence of 0-meson produc- 
tion in proton-nucleus collisions has been carried out. We dis- 
cuss the application of an eikonal formalism, corrections due to 
multiple scattering and the extention to large angle production 
processes. Furthermore, the yl-dependence due to two-step pro- 
duction mechanisms and multi-step processes are investigated 
in detail. We provide all formulas frequently used in the analy- 
sis of nuclear data and study their compatibility and conditions 
of applicability. 

The 0-meson distortion cross section was evaluated from 
the available nuclear data. It was found that different measure- 
ments result in different values of 0$ ranging from to 14 mb. 
Unfortunately, at present the uncertainties of the experimental 
results are too large to draw definite conclusions. 

We also discuss an interpretation of the 0-meson distor- 
tion in nuclear matter and give the relation between various 
frequently used variables, such as in-medium width, distortion 
cross section and scattering length. Furthermore, we show the 
estimates for in the vacuum obtained by VDM, the additive 
quark model, QCD sum rules and other theoretical frameworks 
available. Moreover, we collect predictions for the 0-meson 
modification in matter. While most of the estimates are in rea- 
sonable agreement with <70~1O mb at 0-meson energies above 
3 GeV there are very large uncertainties at lower energies. 

To resolve the unsatisfactory current situation, we propose 
to study the A-dependence of 0-meson production from p+A 
collisions at COSY energies. We provide detailed calculations 
of the ratios of 0-meson production cross sections from dif- 
ferent nuclear targets. Our results can directly be used for the 
evaluation of the from such measurements. 
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